63. 


NOTE ON THE THEORY OF SIMULTANEOUS LINEAR DIF- 
FERENTIAL OR DIFFERENCE EQUATIONS WITH CONSTANT 
COEFFICIENTS. 


[American Journal of Mathematics, Iv. (1881), pp. 321—326.] 


Tus theory is virtually the same for differential as for finite-difference 
equations. The mere verbal part of the exposition „being somewhat easier 
for the former of the two, I shall prefer in the first instance to deal with 
them, although the applications are more interesting when made to bear on 
the latter. Simple to the last degree as are the method of solution and the 
nature of the result, I do not find the one or the other set out, or even 
indicated, except in the most perfunctory manner, in the ordinary text-books. 
This brief notice, designed for the junior readers of the Journal, is intended 
to supply the lacuna. 


Let uj; denote a linear function, with constant coefficients, of wp and of 
its first e; derivatives in respect to t. 


Let Uy,1 + U2 + eee + wù, = 0, 


Us,1 + Uga Fo + Us; = 0, 


Coo eee eee eee rere ee eeeeeeeeeee 


Wia + Uia + 0. F Uii = 0, 
be the system of differential equations proposed for integration. 
Call CRUT LY PRUG 


The process of arriving at the reducing equation for any one of the 
variables is after the manner of the dialytic method of elimination, namely : 
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Along with the first equation take each of its (o — ¢,)th derivatives, with 
the second equation each of its (s — eth derivatives,... and with the ith 
equation each of its (o — ¢;)th derivatives, 


There will thus come into existence (o + 1) i — (e, + e+... +), that is, 
i (e +1)-— c equations between the i(¢+1) quantities 


w, ÒO, eee Àw, 
w2, Òt Os, tee ÒT w2, 

r ° 
wi, twi, .. Ò: Wi. 


If we omit those which appear in any one of the lines above written, there 
will remain (øe + 1)(¢—1) cr i(@ +1)—a@—1 which might be eliminated 
between the i (e + 1) — ø equations, and there would thus result an equation 
between the quantities contained in the omitted line. This elimination, it 
will presently be seen, there is no occasion to perform; the noticeable alge- 
braical fact about it is, that supposing it were performed, the form of the 
equation resulting between wp, Ôt@p, ... &:7@, is invariable, whichever of the 
numbers 1, 2, 3,... 7 be the value assigned to k. 


Let the order of the highest derivative of each œ be reduced by one unit 
below the highest order previously taken, then there will be ic =ø or (i — 1) o 
equations connecting the io quantities 


w, 5:1, 0 tle 5:7 @,, 


W2; Òt w2, eee 5:°1 ws, 


Wi, Òtwi, ose 6:7 @;, 


and accordingly, if we omit the ø quantities which appear in any one (say 
the first) of the above lines, the remaining (¢—1)o¢ quantities may each of 
them be expressed as linear functions of œw, and its (¢ —1) derivatives: but 
the elimination previously indicated would lead to a homogeneous linear 
equation between œw, and its o derivatives, and if in that, each argument 
òw, be replaced by hì and Ay, Az, ... Xe be the o roots of the algebraical 
equation so formed, it follows from the ordinary theory for a single equation 
that œ, (provided the given equations, and consequently the resulting ones, 
be left in their general form) will be of the form 


Cet > a 1C elt T oo at 10g e"o, 


and consequently by virtue of the previous remark @,, œ, ... œp will be of 
the same form as œ, (but, of course, with different coefficients), that is to say, 
the ø roots h, ha, ... ho are the same for the equation in cp as for the equation 
in c, so that the coefficients in the equation between w+ and its ø derivatives 
are, as premised, independent of the value of k. 
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Finally, to determine the equation whose roots are h, ho,... he, let Ce", 
one of the terms in the general value, be taken as a particular value of @,, 
which with corresponding values of the other œs will serve to satisfy the given 
equations; @,, @3,...@; being each of them linear functions of w and deriva- 
tives of w, must be of the forms *Ce”, *Ce*t, ... ‘Ce, so that œ, @,...@; and 
the derivatives of each of them will contain the common factor e*t, and by 
substitution in the original equations we shall obtain a system of simultaneous 
algebraical equations leading to the equation 


Ri; Roy rudy | 
Rar» Tiara ae. LR 


Rio Rie. Rii 
where in general Rp, is what up q becomes on writing h” in place of ô wq. 


The above determinant of the'ith order will be of degree e, + es +... + €i, 
that is, of the degree o (for the general case) in h, and the roots of the equation 
will give the o values A, hy, ... he. 


It follows, therefore, that the result of the hypothetical elimination in the 
first instance referred to will be a linear function of §,%@;, ÒT Op, ... twk, Ok 
of which the coefficients will be identical with the coefficients of h7, he~, ... h, 1 
in the above determinant. Hence no matter now what special values may be 
attributed to the coefficients of the given equations, the result last obtained 
remains of universal validity—without excepting those cases in which the 
result of the hypothetical elimination would be such that the corresponding 
algebraical equation possess equal roots, although in those cases the form 
assumed in the course of the argument for the value of w, (namely, a linear 
function of exponentials) ceases to hold good. Neither for the same reason 
need any exception be made for those cases where the number of terms in 
the equation to wz falls below o on account of one or more of the leading co- 
efficients in the result of the hypothetical elimination becoming zero: the 
degree to which h rises in the determinant will be in all cases the right degree, 
whether it reaches the extreme possible limit o or falls below it. 


The result obtained may be briefly summarized as follows. 
af ($181) @ + (p8) y + -+ + (P81) 2 = 0, 
(Wide) @ + (Yò) y+... + (pid) 2 = 0, 


ee ee ee ee) 


(w) @ + (wt) Y +.. + (wð) 2 =0, 
(each $, W, ..., w standing for a rational-integral functional form) then will 
(Rò) «= 0, (Rò) y= 0, eee (R8,) z = 0, 
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where R (ò) is the resultant in respect to æ, y,...z of what the above 
equations become when © is treated as an ordinary algebraical quantity ; 
under which form the proposition (by virtue of Euler’s method of multipliers) 
becomes so nearly intuitive as to abrogate all necessity for any other 
demonstration *, . 


To pass to the parallel and more important theory in finite differences, it 
is only necessary to interpret uj p to signify a linear function, with constant 
coefficients, of (wr, (@r)lt+ «++ (@x)t+e,> where ¢ is the integer independent 
variable, (say (wp) and its e; difference-augmentatives), and instead of taking 
the differential derivatives of any one of the given equations, to take the 
corresponding difference-augmentatives. Then by precisely the same reason- 
ing as before we shall have 


Oto + Botrert... + Lo; =0, 


B, C, ... L being so taken as that h7 + Bh7-!+...+ L shall be the determinant 
represented by the same form of matrix expressed by R’s as before, but 
where R, , is obtained from u,,¢ by writing h° in lieu of any argument œ: + 0 
which occurs in it. 


The simplest example that can be given is where i= 2, ¢,=e,=1, 
U1 = — Nei + an, Ur, = bhr, 
Wa, 1 = Ct (i = — Orri +d; ; 
this was the case which occurred in the article on the extension of 
Tchebycheff’s theorem, in the last number of the Journal [p. 530, above], 
leading to the equation 
a—Xr b 
c d= 


and to expressions for 7;, 6; as linear functions of Mé, Ax. 


= 0, 


It may also be remarked that this same case gives an instantaneous 
solution of the problem proposed and successfully treated by Babbage in his 
Calculus of Functions, more than half a century ago, and since revived in 
connection with the theory of substitutions (Serret, Alg. Sup. 4 ed., tom. 2, 


pp. 256—262). The problem is to find ġe = pS Aix + Qi 


Dt a i ART 
Bu th so that ‘a, say Fh 
shall equal æ for a given value of 2, 

* I regret that this simple reflection did not present itself to my mind before the preceding 
investigation, the necessity for which it does away with, had been set up in print. It of course 
applies equally well to the analogous proposition for finite-difference equations (u,, vi, ... being 
substituted for x, y ..., and 1+A for ô). This last named proposition, limited to the case of 
equations of the first order, is the foundation-stone of my new theory of Matrices regarded as 
Quantities, that is, as subject to every kind of functional operation which ordinary: arithmetical 
or algebraical quantities are or can be subject to: but though so important and so easily 
established, I know not where it can be found explicitly stated. 
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To find in general ¢‘z it is only necessary to solve the difference equations 
Ui = AU + Aj, 
V; = Buz. + bti, 
and then w, v; will, if u=1, v=0, coincide with a;, Bi, and if m=0, »=1 
with a;, bi. 
Thus calling pı, pa the two roots of 
| —pta a | wn 
B prop 
a; will be of the form Ọ (pt — pŻ) and 8; of the same form except as to O, 
say I'(p;'—p,*). Also a;, b; will be of the forms Opt + Cp, Tipt + Tp, 
where C,+C,=1, T,+T,=1, and the required condition will be fulfilled, 
provided only that p,'=p,', or say 


À AN 
pi= K (cos $T + y(— 1) sin“*) 

A W 
p= K (cos F — y(— 1) sin m 


that is, if (a + b} — 4 (ab — aß) (cos a) =1, A having any integer value 
(which without loss of generality may be taken inferior to 7) except zero*. 
If à= 0, the two roots of the equation in p become equal and the form of 
the solution changes into 
= (C,+ Ori) p', 14 = (0 + O10) pt 
When w= 1 and v,=0 then u =a, » = £, 
C,=1, G'=0, 0 =5-1, joa, 


p 
and when w= 0, v=1l1, w=a, n= b, 


Qat, akii y C= 0, C,=1, 
p p 


and ¢'a = te a yee = , which cannot be pence for any value of t, 


and when 7 = œ becomes 


(a—p)@+a_a-—p_ a a-b 2a 
Bet+b—p B b-p’ 28 a—b’ 


so that tæ in this case continually converges to a constant limit. 


that is, = 


I may add that ¢’# converges to a constant limit not merely when the 
roots pı, pa of 
a-—p a 
| B a P 
* There will thus be (i — 1) values of \ which will each give a distinct admissible solution of 
the problem of periodicity, but of course only those values of \ which are relatively prime to i will 


give primitive solutions. If i=i’5 the effect of making \=)’é will be to make tx =< by virtue o 
its making ix =0. 
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are equal, but whenever they are real. For the general form of ¢‘z, it may 
easily be found, is 

[(p2— @) pr! + (p: — @) pi] £ +.4(p,'— ps") 

B (pi! — po’) @ + [(p2— b) pr’ + (pi — b) p] 
(pi—a)a@—a a-—p, a 
—Be+p,—b 8B b-p 
pı signifies the smaller of the two roots pı, pa; or in other words when 
ax +a 


Bath” is 
, With the understanding that the quantity under 


where 


which if p.>p,; when i=% becomes 


a — b > 2 y(aß), the limiting value to ¢‘x, when æ represents 
(a — b) +v {(a — b} — 4ap) 
2 


the radical sign is to be taken positive. 
So, if 
Zip | Yip | Sign = UL; + by: + C2; : Wat byi + o'zi : aw; + by; + ezi, 
when all the roots of the determinant 
| a-Hxr b c 
a b- c 
| a” p shi a X 
are real, the point a;, Yı Zi, as i increases, will be found to approach 


indefinitely near to a fixed straight line; and if all the roots are equal, to a 
fixed point. 


The condition of the system of ratios æ; :y;:z2; being periodic and 
having a period m is tantamount to the condition that the mth power of 
the matrix 


rae | ae 
a Pr 
o E ae hg 
shall be the matrix 
T i 
ORRIA 
Di Oy ks 


The complete solution of this problem, and of the more general one of 
extracting the mth root of any unit-matrix (that is, a matrix in which each 
element in the principal diagonal is unity, and the rest zero), which constitutes 
the ultimate generalization of Babbage’s problem and is soluble by the same 
method, will probably appear in a memoir on matrices, in the forthcoming 
number of the Journal. 


In general, for a matrix of the order w, the number of mth roots is m® and 
each of them is perfectly determinate. But when the matrix is a unit-matrix 
or a zero-matrix (the latter meaning one in which every element is zero) 
there are distinct genera and species of such roots, and every species contains 
its own appropriate number of arbitrary constants. 
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